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Abstract 


A  Note  on  Programming  and  the  Single  Machine 
Lotsize  Scheduling  Problem 

by 

Thom  J.  Hodgson 
Henry  L..  W.  Nuttle 


Past  research  on  the  lotsize  scheduling  problem  has  relied  heavily 
on  the  following  assumptions: 

1.  Each  product  is  to  be  produced  on  a  regular  invariant  cycle. 

\ 

2.  The  inventory  level  at  the  onset  of  each  production  run  is  zero. 

In  this  paper  these  two  restrictions  have  been  dropped.  It  is  shown  that  the 
inventory  cost  function  can  be  expressed  as  a  linear  function.  This  allows  the 


use  of  linear  programming  methods  to  determine  the  optimal  lengths  for  a  given 
sequence  of  production  runs.  'JL-  ■ 
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The  lotsize  scheduling  problem  is  basically  that  of  finding  a  minimum  cost 
production  schedule  on  a  single  machine,  for  N  different  products,  each  with  known 
demand  rate,  production  rate,  and  setup  time.  As  its  name  implies,  the  problem  has 
elements  of  both  inventory  and  scheduling.  The  objective  function,  the  sum  of  setup 
and  inventory  costs,  is  derived  from  inventory  theory;  and  the  feasibility  require¬ 
ment  associated  with  machine  availability  is  primarily  a  scheduling  problem. 

Variations  of  the  problem  have  been  the  subject  of  numerous  articles  in  the 
literature.  A  comprehensive  bibliography  through  1976  is  given  in  [4].  More  recent 
references  include  [ 1 ] — [ 3 ] , [ 5 ] — [ 7 ] ,  [9],  and  [11]. 

The  most  frequently  considered  variation  is  the  one  in  which  the  objective  is 
to  determine  a  cyclic,  repeatable  schedule  which  minimizes  setup  plus  inventory 
costs  per  unit  time.  In  this  context  one  is  faced  with  determining  the  sequence  in 
which  products  are  to  be  processed  and  the  length  of  each  production  run. 

In  this  paper  we  focus  on  determining  the  run  lengths,  given  the  sequence.  In 
particular,  we  show  that  when  all  runs  for  a  given  product  must  be  of  equal  length, 
but  production  may  start  before  inventory  is  exhausted,  then  the  optimal  cycle  and 
run  lengths  may  be  determined  by  linear  programming  methods.  This  result  can  be 
used  in  conjunction  with  a  sequence  generating  procedure  to  provide  solutions  to  the 
overall  problem.  It  can  also  be  used  to  easily  "tune"  schedules  produced  by  other 
procedures  to  be  more  acceptable  from  the  standpoint  of  implementation. 

Our  model  corresponds  to  Maxwell's  [8]  when  the  "equal  lots  rule"  is  used  and 
idle  time  between  the  production  runs  is  permitted.  Relative  to  the  model  of  Del- 
porate  and  Thomas  [3],  it  requires  equal  lots  for  each  run  of  a  given  product,  but 
allows  non-zero  inventories  at  the  beginning  of  the  runs. 
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Consider  a  set  of  N  products  to  be  produced  in  lots  according  to  a  repeating 
cyclic  schedule  on  a  single  machine  which  can  process  one  product  at  a  time.  Let 
d^,  pk,  and  Sk  denote  demand  rate,  production  rate,  and  set  up  time,  respectively, 
for  product  k. 

Suppose  that  we  are  given  a  sequence  J- {k^ , . . ,kj , . . ,k^  }  of  production  runs  to 

be  performed  in  a  cycle  of  length  T  (to  be  determined)  tine  units. 

k j e { 1 , 2 , . . . , iJ  }  denotes  the  product  to  be  processed  during  run  j  and  n^  denotes  the 

cardinality  of  tbg  set  {  j  e[  1 , 2 , . . . ,  L]  |  k^k),  i.e.,  the  number  of  runs  of  product  k. 
N  J 

Thus  L  =  £  n^.  We  require  all  production  runs  for  product  k  to  yield  d^T/n^  units 

k=l 

which  implies  a  run  length  of  d^T/p^n^  time  units. 

In  this  context,  we  wish  to  determine,  the  cycle  length,  T,  and  the  production 
run  start  times  (after  setup),  Xj,  0  ^  Xj  ^  T,  j  =  1,  ...,  L,  so  as  to  minimize 
the  cost  (setup  plus  inventory  carrying)  per  unit  time,  subject  to  the  restriction 
that  all  demand  be  satisfied  without  backordering. 

If  Wj  denotes  the  amount  of  idle  time  following  production  run  j,  then  start 
times,  idle  times ,  and  cycle  length  are  related  by 

dk  T 
j 

Xj  + - +  wj  +  sk  =  X j+1  *  j  =  1,  ...,  L  -  1  (1) 

Pk  "k  j+1 

j  j 

dk  T 
L 

XL  + - +  \  =  T 

Pk  n:< 

L  L 

where  X^  =  . 

1 

Let  Zj  be  the  index  of  the  next  production  run  of  product  kj  after  production 

run  j.  If  1 j  is  the  inventory  level  of  product  kj  at  the  beginning  of  run  j,  then 

U  +  d.T/n,  -  d,(X  -  X.)  =  I  ,  j  =  1,  2,  ...  (3) 

j  k.  ^  ^  J  zj 


3 


along  with  I j  0,  for  all  j,  assures  that  no  shortages  will  occur.  Note  that  equa¬ 
tion  (3)  is  not  neeeded  when  production  run  j  is  the  n^th  production  run  for  product 


j 

kj,  since  only  nk  -  1  equations  are  necessary  to  fully  define  the  relationships  be- 

j 

tween  the  nk  inventory  levels. 


j 


If 


Our  objective 
inventory  levels 


is  to  minimize  the  total  setup  plus  carrying  cost  per  unit  time. 
,  I j ,  were  required  to  be  zero,  then  this  cost  would  be  given  by 


N 

Z  (nkAk/T  +  hck(l-dk/pk)dkT/2nk)  (4) 

k=! 


where  A,  and  c.  are  the  setup  and  piece  costs,  respectively,  for  product  k,  and  h  is 

the  inventory  carrying  charge  .  The  associated  inventory  pattern  for  an  arbitrary 

product  k  (nk  =  3)  is  shown  in  Figure  la.  When  the  Ij's  are  allowed  to  be  positive 

(implying  that  production  of  a  product  may  begin  before  the  current  supply  is 

exhausted)  inventory  patterns  such  as  that  illustrated  in  Figure  lb  (where  we  use 

the  notation  Ij  to  represent  the  level  for  product  k  at  the  beginning  of  its  ith 
i 

run)  are  then  possible.  Such  patterns  result  in  an  increased  carrying  cost  (versus 

I  ;  =  0)  which  is  proportional  to  the  area  (A+B)  in  Figure  lb  (or  its  equivalent  for 
i 

other  patterns).  Since  the  production  schedule  will  be  repeated  every  T  time  units, 
the  increase  in  carrying  cost  for  product  k  per  unit  time  may  be  calculated  as 
hck(A+E)/T  (or  its  equivalent  for  other  patterns).  As  illustrated  in  Figure  lc, 
this  increase  is  easily  shown  to  be 


hck  1  Xj  /nk 
i  =  l  i 


(5) 


Note  that  the  increase  is  linearly  dependent  on  the  minimum  inventory  levels. 
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Based  on  equation  (5),  the  cost  expression  ((4))  may  be  rewritten  as 

N  N  L 

l  nkVT  +  ht  1  ck(1-dk/Pk)dkT/2nk  +  E  ck  ri/nk  1  <6> 

k=l  k=l  j=l  J  J 

Minimizing  (6)  subject  to  (1)  -  (3)  plus  nonnegativity  restrictions  on  the  X^'s, 

Ij's,  Wj's  and  T  will  yield  an  optimal  schedule. 

This  problem  may  be  solved  by  linear  programming  methods.  Note  that  the 
problem  may  be  restated  as,  choose  T  >  0  to 

min  Z(T)  =  A(T)+F(T)  (7) 

N 

where  A(T)  =  £  n.  A,  /T  and 

k-1  ^  k 

N  L 

F(T)  =  min  h  [  l  ck( 1-dk/Pk)dkT/2nk  +  Z  ck  Ii/nk  1 
k=l  "  j  =  l  j  j 

subject  to  ( l )  —  ( 3 )  plus  nonnegativity,  which  is  a  linear  program. 

For  T>0,  the  hyperbolic  setup  cost  function,  A(T),  is  convex,  while  F(T)  is 
piecewise  linear  and  convex  [10].  Thus  Z(T)  is  convex  for  T>0. 

Let  T°  be  the  point  where  F(T)  attains  its  minimum.  For  T<T°  both  A(T)  and 
F(T),  and  thus  Z(T),  are  decreasing.  Thus  the  optimum  T  will  be  at  least  as  large 
a:;  T° ,  so  that  we  can  restrict  our  attention  to  T>T° .  For  this  range  F(T)  can  be 
obtained  by  adding  the  constraint 

T  >  T°  (8) 

to  the  linear  program  and  applying  sensitivity  analysis  to  T°.  For  T  >  T° ,  F(T)  is 
piecewise  linear  and  non  decreasing  with  break  points  corresponding  to  optimal  basis 
changes  in  the  Li’.  Thus,  the  minimum  of  (7)  will  occur  either  at  one  of  the  break- 
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points  of  F(T)  or  in  between  breakpoints.  To  determine  if  the  minimum  occurs  be¬ 
tween  two  breakpoints,  the  following  procedure  can  be  used.  The  expression  for  the 
variable  elements  of  the  total  cost  function  in  a  given  segment  is  just 


£  n.  A,/T  +  ST, 
k=l 


where  S  is  the  sensitivity  to  T  (dual  variable  for  (8))  of  the  linear  program  objec¬ 
tive  function  in  the  segment.  The  minimum  of  (9)  occurs  at 


Z  nkAk 


If  T'  falls  within  the  boundaries  of  the  segment,  then  the  minimum  of  Z(T)  occurs  at 
T'.  Since  S  is  non  decreasing  with  T° ,  if  T'  falls  below  the  lower  end  of  the 
segment,  then  the  optimal  T  occurs  at  the  lower  end  of  the  segment.  If  T'  falls 
above  the  upper  end  of  the  segment,  the  next  segment  must  be  investigated. 

If  we  ignore  the  requirement  that  only  one  product  be  processed  at  a  time,  the 
least  cost  value  of  T  can  be  determined  as 


2  E  nkAk 


h  L  ckdk(1~Vpk>/nk 

l  k=1 

When  this  requirement  is  considered,  it  is  easy  to  see  that  the  optimal  value  of  T 
will  be  no  less  than  T* .  Thus  we  may  replace  constraint  (8)  with 


T>raax  (T° ,  T*) 


Finally,  we  should  note  that,  from  the  standpoint  of  implementation,  the 
production  quantities  associated  with  the  optimal  value  of  T  may  be  impractical. 
For  example,  if  T  turns  out  to  be  19.335  workdays,  this  may  call  for  producing  2 
lots  of  a  product  each  of  which  covers  the  demand  for  9.6675  workdays.  As  an 
alternative,  albeit  at  some  cost  relative  to  optimality,  we  may  simply  specify  a  T" 

(  for  example,  20  workdays)  and  solve  the  LP  with  T  =  T” .  In  this  context,  the  op¬ 
timal  T  can  be  used  to  determine  the  "ballpark"  for  the  choice  of  T" . 
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